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The flow of fluids in channels, pipes or ducts, as in any other wall-bounded flow (like water along 
the hulls of ships or air on airplanes) is hindered by a drag, which increases many-folds when 
the fluid flow turns from laminar to turbulent. A major technological problem is how to reduce 
this drag in order to minimize the expense of transporting fluids like oil in pipelines, or to move 
ships in the ocean. It was discovered in the mid-twentieth century that minute concentrations of 
polymers can reduce the drag in turbulent flows by up to 80%. While experimental knowledge 
had accumulated over the years, the fundamental theory of drag reduction by polymers remained 
elusive for a long time, with arguments raging whether this is a "skin" or a "bulk" effect. In 
this colloquium review we first summarize the phenomenology of drag reduction by polymers, 
stressing both its universal and non-universal aspects, and then proceed to review a recent theory 
that provides a quantitative explanation of all the known phenomenology. We treat both flexible 
and rod-like polymers, explaining the existence of universal properties like the Maximum Drag 
Reduction (MDR) asymptote, as well as non-universal cross-over phenomena that depend on the 
Reynolds number, on the nature of the polymer and on its concentration. Finally we also discuss 
other agents for drag reduction with a stress on the important example of bubbles. 
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I. INTRODUCTION 

The fact that minute concentrations of flexible poly- 
mers can reduce the drag in turbulent flows in straight 
tubes was first d i scover ed by B.A. Toms and published 
in 1949 (jToms I . ll949T ). This is an important phe- 
nomenon, utilized in a number of technological applica- 
tions including the Trans-Alaska Pipeline System. By 
1995 there were about 2500 papers on the subject, and 
by now t here are many more. Earlier reviews include 
those by (Ide-Gennest 1 19901: iHovtl. Il972t iLandhall. 11975 ; 
Lumlevl. 119691: iMcCombT 1990; Src enivasan and W hite. 
20001 : IVirkl " I1975T ). and others. In this introductory chap- 
ter we present some essential facts known about the phe- 
nomenon, to set up the theoretical discussions that will 
form the bulk of this review. 

It should be stated right away that the phenomenon 
of drag, which is distinguished from viscous dissipation, 
should be discussed in the context of wall-bounded flows. 
In hom ogeneous isot ropic turbulence there exists dissi- 
pation (jLambl . I1879I ). which at every point of the flow 
equals i^o | Vi7(r, t)\ 2 where Vq is the kinematic viscosity 
and U(r, t) is the velocity field as a function of position 
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FIG. 1 The channel geometry 



and time. The existence of a wall breaks homogeneity, 
and together with the boundary condition U — on the 
wall it sets a momentum flux from the bulk to the wall. 
This momentum flux is responsible for the drag, since not 
all the work done to push the fluid can translate into mo- 
mentum in the stream-wise direction. In stationary con- 
dition all the momentum produced by the p ressure head 
must flow to the wall. Understanding this (jL'vov et all 
l2004j ) is the first step in deciphering the riddle of the 
phenomenon of drag reduction by additives, polymers or 
others, since usually these agents tend to increase the 
viscosity. It appears therefore counter-intuitive that they 
would do any good, unless one understands that the main 
reason for drag-reduction when polymers or other drag- 
reducing agents are added to a Newtonian fluid is caused 
by reducing the momentum flux to the wall. The rest of 
this review elaborates on this point and makes it quanti- 
tative. 

Contrary to the preference of the engineering commu- 
nity, in which arguments can rage about what is the 
mechanism of drag reduction, and whether these are hair- 
pin vortices or other 'things' that do the trick, we will 
adhere to the parlance of the physics community where 
a theory is tested by its quantitative predictiveness. We 
thus base our considerations on the analysis of model 
equations and their consequences. The criterion for va- 
lidity will be our ability to describe and understand in a 
quantitative fashion all the observed phenomena of drag 
reduction, both universal and non-universal, and the abil- 
ity to predict the results of yet unperformed experiments. 
The reader will judge for himself the extent of success of 
this approach. 



A. Universality of Newtonian mean velocity profile 

For concreteness we will focus on channel flows, but 
most of our observations will be equally applicable to 
other wall bounded flows. The channel geometry is 
sketched in Fig. [1] The mean flow is in the x direc- 
tion, between two parallel plates displaced by a distance 
2L. The distance from the lower wall is denoted by y, 
and the span-wise direction by z. One takes the length of 
the channel in the x direction to be much larger than the 
distance between the side walls in z, and the latter much 
larger than L. In such a geometry the mean velocity 



V = (U(r, t)) is (to a high approximation) independent 
of either x or z, but only a function of the distance from 
the wall, V = V(y). When a Newtonian fluid flows at 
large large Reynolds number (cf. Eq. in such a 

channel, it exhibits in the near-wall region a universal 
mean velocity profile. Here we use the word "universal" 
in the sense that any Newtonian fluid flowing in the vicin- 
ity of a smooth surface will have the same mean velocity 
profile when plotted in the right coordinates. The uni- 
versality is best displayed in dimeniso nless coordinates, 
known also as "wall units" (|PopeLl2000T l; First, for incom- 
pressible fluids we can take the density as unity, p = 1. 
Then we define the Reynolds number Tie, the normal- 
ized distance from the wall y + and the normalized mean 
velocity V + (y + ) as follows: 



Tie = L^p'L/ Vo , y + = yTie/L , V + = V/ 



(1.1) 



Here p' is the fixed pressure gradient p' = ~dp/dx. The 
universal p rofile is shown in Fig. [21 with green circles [sim- 
ulations of lDe-Angelis et a l. (2003)1 open cir cles (experi- 
ments ) and a black continuous line [theory by lL'vov et all 
(2004)]. This profile has two distinct parts. For y + < 6 
one observes the "viscous sub-layer" where 



V + (y+)=y 



y 



< 6 



(1.2) 



(and see Subsect. ll.3l for a derivation), whereas for y + > 
30 one sees the celebrated universal von Karman "log-law 
of the wall" which is written in wall units as 



V + (y + ) = n- 1 \ny+ + B K , for y+ > 30 



(1.3) 



The law (|1.3p is universal, independent of the nature of 
the Newtonian fluid; it had been a shortcoming of the 
theory of wall-bounded turbulence that the von Karman 
constant k k w 0.436 and the intercept B K w 6.13 
had been only know n from experiments and sim u lation s 
(|Monin and Yagloml . Il979t IZagarola and Smitsl. I1997T). 
Some recent progress on this was made by lLo et al.l 
(120051) . 

Having observed the universal profile for the mean 
velocity, it is easy to see that any theory that seeks 
to understand dr ag reduction by a change in viscosity 
(Ide-Gennesl . Il990h is bound to fail, since the universal 
profile is written in reduced coordinates, and a change in 
the viscosity can result only in a re-parameterization of 
the profile. Indeed, drag reduction must mean a change 
in the universal profile, such that the velocity V + in re- 
duced coordinates exceeds the velocity V + predicted by 
Eq. (HI. 



B. Drag reduction phenomenology 

Here we detai l some of the prominent features of drag 
reduction ( Virkl . Il975h . all of which must be explained 
by a consistent theory. An immediate riddle that comes 
to mind is the following: the polymers are molecular in 
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FIG. 2 Mean normalized velocity profiles as a function of the 
normalized distance from the wall. T he data points (gr e en cir - 
cles) from numerical simulations of iDe-Angelis et al.l (I2003T) 
and the experimental points (open circles dWarholic et all 
1 19991 ) refer to Newtonian flows. The s olid li ne is a theoret- 
ical formula deve loped by iL'vov et al. | fill). The red data 
points (squares) l|Virkl I1975T ) represent the Maximum Drag 
Reduction (MDR) asymptote. The dashed red curve repre- 
sents t he log-law (11 . 5 p wh ich was derived from first princi- 
ples b y lBenzi et al] (|2005ah . The blue fi l led tr iangles (|Rollinl. 
I1972T ) and green open triangles l|Ruddl 1 1969T > represent the 
cross-over, for intermediate concentrations of the polymer, 
from the MDR asymptote to the Newtonian plug. 



scale. Turbulence is characterized by the outer scale L 
where energy is injected into the system, and by the Kol- 
mogorov viscous scale r\ below which viscous dissipation 
dominates over inertial terms, and the velocity field be- 
comes essentially smooth. For all realistic flows the poly- 
mer size is much smaller than this viscous scale (and 
cf. Subsec. IIII.AI for some actual numbers for a typical 
polymer). How is it then that the polymers can inter- 
act at all with the turbulen t degre es of freedom? This 
riddle was solve by iLumlevI (jl969| ) who argued that it 
is the polymer relaxation time r, the time that charac- 
terizes the relaxation of a stretched polymer back to its 
coiled equilibrium state, which is comparable to a typi- 
cal eddy turn-over time in the turbulent cascade. This 
matching of time scales allows an efficient interaction be- 
tween the turbulent fluctuations and the polymer degrees 
of freedom. With the typical shear rate S(y) one forms 
a dimensionless "Deborah number" 



Ve{y) = rS(y) 



(1.4) 



When Ve exceeds the order of unity, the polymers begin 
to interact with the turbulent flow by stretching and tak- 
ing energy from the turbulent fluctuations (cf. Subsec. 
IIII.AI for a derivation of this). We will see below that this 
mechanism of Lumley is corroborated by all the available 



data. What we will have to explain is how is it that as 
the polymers stretch, a process that must increase the 
viscosity, nevertheless the drag reduces. 



1. The universal Maximum Drag Reduction asymptote 

One of t he most significant experimental findings 
(|Virkl . I1975D concerning turbulent drag reduction by 
polymers is that in wall-bounded turbulence (like chan- 
nel and pipe flows) the velocity profile (with polymers 
added to the Newtonian fluid) is bounded between the 
von Karman's log-law (|1.3|) and another log-law which de- 
scribes the maximal possible velocity profile (Maximum 
Drag Reduction, MDR), 



(1.5) 



where k^ 1 w 11.7 and B v w —17. This law, which 
was discovered experimentally bv lVirkl ( 19751 ) (and hence 
the notation k v ), is also claimed to be universal, inde- 
pendent of the Newtonian fluid and the nature of the 
poly mer additive, incl uding flexible and rod-like poly- 
mers ()Virk et allll997D . This log-law, like von Karman's 
log-law, contains t wo phenomenological parameters. In 
( L'vov et al. , 2004) it was shown that in fact this law 
contains only one parameter, and can be written in the 
form: 

V+(y+) = K- 1 \n(eK v y+) for y+ > 12 , (1.6) 



where e is the basis of the natural logarithm. The deep 
reason for this simplification will be explained below. 
For sufficiently high values of Tie, and sufficiently high 
concentration of the polymer c p , and length of polymer 
(number of monomers N p ), the velocity profile in a chan- 
nel is expected to follow the law (|1.6|) . Needless to say, 
the first role of a theory of drag reduction is to provide an 
explanation for the MDR law and for its universality. We 
will explain below that the deep reason for the universal- 
ity of the MDR is that it is a marginal state between a 
turbulent and a laminar regime of wall bounded-flows. 
In this marginal state turbulent fluctuations almost do 
not contribute to the momentum and energy balance, 
and the only role of turbulence is to extend the polymers 
in a proper way. This explanation can be found in Sect. 
IIII1 including an a- priori calculation of the parameter k v . 
For finite Tie, finite concentration c p , and finite number 
of monomers N p , one expects cross-overs that are non- 
universal; in particular such cross-overs depend on the 
nature of the polymer, whether it is flexible or rod-like. 



2. Cross-overs with flexible polymers 

When the drag reducing agent is a flexible polymer, 
but the concentration c p of the polymer is not sufficiently 
large, the mean velocity profiles exhibit a cross-over back 
to a log-law which is parallel to the law (|1.3[) . but with a 
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larger mean velocity (i.e. with a larger value of the inter- 
cept B K ), see Fig. [2] The region of this log-law is known 
as the "Newtonian plug" . The position of the cross-overs 
are not universal in the sense that they depend on the 
nature of the polymers and the flow conditions. The sce- 
nario is that the mean velocity profile follows the MDR 
up to a certain point after which it crosses back to the 
Newtonian plug. The layer of y + values between the vis- 
cous layer and the Newtonian plug is referred to as the 
"elastic layer" . A theory for these cross-overs is provided 
below, cf. Sect. El 

Another interesting experimental p iece of info r matio n 
about cross-overs was provided by IChoi et al (2002). 
Here turbulence was produced in a counter-rotating disks 
apparatus, with A-DNA molecules used to reduce the 
drag. The Reynolds number was relatively high (the re- 
sults below pertain to Re « 1.2 x 10 6 ) and the initial 
concentrations c p of DNA were relatively low [results em- 
ployed below pertain to c p — 2.70 and c p — 1.35 weight 
parts per million (wppm)]. During the experiment DNA 
degrades; fortunately the degradation is very predictable: 
double stranded molecules with 48 502 base pairs (bp) in 
size degrade to double stranded molecules with 23 100 
bp. Thus invariably the length N p reduces by a factor of 
approximately 2, and the concentration c p increases by 
a factor of 2. The experiment followed the drag reduc- 
tion efficacy measured in terms of the percentage drag 
reduction defined by 
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FIG. 3 %DR in a counter-rotating disks experiment with A- 
DNA as the drag reducing polymer. Note that the %DR is 
proportional to c p . When the length N p reduces by a factor 
of 2 and, simultaneously, c p increases by factor of 2, the %DR 
reduces by a factor of 4. 



%DR = Tn t Tv x 100 . 



(1.7) 



where T N and T v are the torques needed to maintain 
the disk to rotate at a particular Reynolds number with- 
out and with polymers, respectively. The main exper- 
imental results which are of interest to us are summa- 
rized in Fig. [3l We see from the experiment that both 
initially (with un-degraded DNA) and finally (with de- 
graded DNA) the %DR is proportional to c p . Upon de- 
grading, which amounts to decreasing the length N p by a 
factor of approximately 2 and, simultaneously increasing 
c p by factor of 2, %DR decreases by a factor 4. Expla- 
nations of these findings can be found in Subsect. IV.BI 



3. Cross-overs with rod-like polymers 

For the purposes of this review a flexible polymer is 
a polymer that is coiled at equilibrium or in a flow of 
low Reynolds number, and it undergoes a coil-stretch 
transition at some value of the Reynolds number, (see 
below for details). A rod- like polymer is stretched a- 
priori, having roughly the same linear extent at any value 
of the Reynolds number. When the polymers are rod- 
like and the concentration c p is not sufficiently high, the 
cross-over scenario is different. The dat a in Fig. [4] in- 
clude both flexible and rod-like polymers (|Escudier et all 
1999). It again indicates that for large values of Tie the 
mean velocity profile with flexible polymers [polyacry- 



lamide (PAA)] follows the MDR until a point of cross- 
over back to the "Newtonian plug" , where it becomes 
roughly parallel to von Karman's log-law. Increasing 
the concentration results in following the MDR further 
until a higher cross-over point is attained back to the 
Newtonian plug. On the other hand, for rod-like poly- 
mers [sodium carboxymethylcellulose (CMC) and sodium 
carboxymethylcellulose/xanthan gum blend (CMC/XG)] 
the data shown in Fig. [4] indicate a different scenario. 
Contrary to flexible polymers, here, as a function of the 
concentration, one finds mean velocity profiles that inter- 
polate between the two asymptotes (|1.3[) and (|1.5[) . reach- 
ing the MDR only for large concentrations. A similar dif- 
ference in the behavior of flexible and rod-like polymers 
when plotting th e drag as a f unctio n of Reynolds number 
was reported by IVirk et al.l ( 1997? ). Clearly, an explana- 
tion of these differences between the way the MDR is 
attained must be a part of the theory of drag reduction, 
and cf. Sect. IVII We reiterate that these cross-overs 



pertain to the situation in which Tie is large, but c p is 
too small. 

Another major difference between the two classes of 
polymers is found when Tie is too small, since then rod- 
like polymers can cause drag enhancement, whereas flex- 
ible polymer never cause drag enhancement. The lat- 
ter are either neutral or cause drag reduction. The best 
way to see the phenomenon of drag enhancement at low 
Tie with rod-like polymers is to consider the Fanning drag 
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FIG. 5 The drag in Prandtl-Karman coordinates for a low 

concentration solution of NaCl in water mixed w ith the rod- 

10 100 1000 like p olymer PAMH B1120 in a pipe flow, see (Virk et al.l . 

y* I1996T ) for details. The symbols represent the concentrations 

in wppm (weight parts per million) as given in the figure. 



FIG. 4 Typical velocity profiles taken from l|Escudier et all . 
Il999f ). In dashed line we noted the von Karman law (|1.3|l . 
while the MDR (|f .6[l is shown as the continuous black line. In 
all cases the mean velocity follows the same viscous behavior 
for y + < 10. After that the scenario is different for flexible 
and rod-like polymers. The typical behavior for the former 
is presented by the open triangles, which follow the MDR up 
to a cross-over point that depends on the concentration of 
the polymer and on the value of Re. The rod-like behavior is 
exemplified by the solid triangles and the open squares; the 
mean velocity profiles appear to interpolate smoothly between 
the two asymptotes as a function of the concentration of the 
rod-like polymer. 



coefficient denned as 



/-*/(l,v>). 



(1.8) 



where r* is the shear stress at the wall, determined by 
the value S(y = 0) of the shear at y = 0: 



pv S(y = 0) , 



(1.9) 



p and V are the fluid density and the mean fluid through- 
put, respectively. Fig. [5] presents this quantity as a func- 
tion of Tie in the traditional Prandtl-Karman coordi- 
nates 1/V7 vs - T&VJ) f° r which once again the New- 
tonian high Tie log-law is universal, and for which also 
there exists and MDR universal maximum (|VirkL 119751: 
IVirk et al ., 19%). The straight continuous line denoted 
by 'N' presents the Newtonian universal law. Data points 
below this line are indicative of a drag enhancement, i.e., 
an increase in the dissipation due to the addition of the 
polymer. Conversely, data points above the line corre- 
spond to drag reduction, which is always bound by the 
MDR asymptote represented by the dashed line denoted 
by 'M'. This figure shows data for a rod- like polymer 
(a polyelectrolyte in aqueous solution at very low salt 
concentration) and shows how drag enhancement for low 



values of Re crosses over to drag reduc t ion at large values 
of Re (jVirk et all . [l996tl Wagner et al.l . 120031 ). One of the 
results of the theory presented below is that it reproduces 
the phenomena shown in Fig. [5] in a satisfactory manner. 
We reiterate that with flexible polymers the situation is 
very different, and there is no drag enhancement at any 
value of Re. The reason for this distinction will be made 
clear below as well. 



II. SIMPLE THEORY OF THE VON KARMAN LAW 

As an introduction to the derivation of the MDR 
asymptote for fluids laden with polymers we remind the 
reader first how the von Karman log-law (| 1 . 3[) is derived. 
The derivation of the MDR will follow closely similar 
ideas with the modifications due to the polymers taken 
carefully into account. 

Wall- bounded turbulence in Newtonian flui ds is dis- 
cussed (|Monin and Yagloml 1 19791 : iPopel l2000l) by con- 
sidering the fluid velocity U(r) as a sum of its average 
(over time) and a fluctuating part: 



U(r,t) = V(y)x + u(r,t) 



(2.1) 



The objects that enter the Newtonian theory are the 
mean shear S(y), the Reynolds stress W(y) and the ki- 
netic energy density per unit mass K(y) : 

S(y) = dV(y)/dy, W (y) = - (u x u v ) , K (y) = (\u\ 2 )/2 . 

(2.2) 

Note that the Reynolds stress is nothing but the momen- 
tum in the stream-wise direction x transported by the 
fluctuations u y in the direction of the wall; it is the mo- 
mentum flux to the wall. Using the Navier-Stokes equa- 
tions one can calculate this momentum flux P(y) which 
is generated by the pressure head p' = —dp/d x\ at a 
distance y from the wall this flux is (Pop3, l2000l ) 

P(y)=p'(L-y) . (2.3) 
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We will be interested in positions y <C L, so that we can 
approximate the production of momentum by a constant 
p'L. This momentum is then dissipated at a rate v$S{y). 
This give rise to the exact momentum balance equation 



v S{y) + W(y) =p'L, y « L 



(2.4) 



The two terms on the LHS take their main roles at differ- 
ent values of y. For y close to the wall the viscous term 
dominates, predicting that S(y) ~ p'L/vq. This trans- 
lates immediately to Eq. (|1.2p . the mean velocity profile 
close to the wall. 

Away from the wall the Reynolds stress dominates, but 
all that can be learned is that W(y) ^constant, which is 
not enough to predict the velocity profile. We need now 
to invoke a second equation that describes the energy 
balance. 

Directly from the Navier-Stokes equations one can 
compute the rate of turbulent kinetic energy produc- 
tion b y the mean shear, it is W(y)S(y) (see for ex- 
ample \Poni (|2000l) ). and the energy dissipation Ed at 
any point, Ed — i/ (\Vu(r, t)\ 2 ). The energy dissipa- 
tion is estimated d ifferently near the wall and in the bulk 
( L'vov et al.Ll2004j) . Near the wall the velocity is smooth 
and we can estimate the gradient by the distance from 
the wall, and thus Ed ss VQa 2 K(y)/y 2 where a is a dimcn- 
sionlcss coefficient of the order of unity. Further away 
from the wall the flow is turbulent, and we estimate the 
energy flux by K(y)/r(y) where r(y) « y/b-\/K(y) is the 
typical eddy turn-over time at a distance y from the wall, 
and b is another dimensionless coefficient of the order of 
unity. Putting things together yields the energy balance 
equation 



y 



K(y) = W{y)S{y) , 



(2.5) 



where the first term is dominant near the wall and the 
second in the bulk. The interpretation of this equation is 
similar to the momentum balance equation except that 
the latter is exact; the first term on the LHS is the viscous 
dissipation, the second the energy flux to the wall, and 
the RHS is the production. 

To close the problem, the balance equations need to 
be supplemented by a relation between K(y) and W(y). 
Dimcnsionally these objects are the same, and therefore 
in the bulk [where one uses dimensional analysis to derive 
the second term of (|2.5|) ]. we expect that these objects 
must be proporti onal to each oth e r. Indeed, experim ents 
and simulations (Ptasinski ct al 1 120011 IVirkL I1975T) cor- 
roborate th is expectation and one writes ( L'vov et all 
l2004l2005cD 



W{v)=<*K{y) 



(2.6) 



with c N being apparently y-independent outside the vis- 
cous boundary layer. To derive the von Karman log- 
law we now assert that in the bulk the first term on the 
LHS of (|2.5p is negligible, we then use Eq. (|2.6p together 



with the previous conclusion that W(y) =const. to de- 
rive immediately S(y) oc 1/y. Integrating yields the von 
Karman log-law. Note that on the face of it we have 
three phcnomcnological parameters i.e. a, b and c N . In 
fact in the calculation of the mean velocity they appear 
in two combination which can be chosen as k k and B K 
in Eq. (|1.3p . Nevertheless it is not known how to com- 
pute these two parameters a-priori. Significantly, we will 
argue below that the slope of the MDR can be computed 
a-priori, cf. Subsect. IIII. 



III. THE UNIVERSALITY OF THE MDR 

To study the implications of adding small concentra- 
tions of polymers into wall-bounded turbulent flows we 
need a model that describes reliably the modifications 
in the hydrodynamic equations that are induced by the 
polymers. In the present section we are aiming at the uni- 
versal, model independent properties of drag reduction, 
and therefore any reasonable model of polymers interact- 
ing with hydrodynamics which exhibits drag reduction 
should also lead to the universal properties. We begin 
with the case of flexible polymers. 



A. Model equations for flows laden with flexible-polymers 

To give a flavor of the type of polymers most com- 
monly used in experiments and technological applications 
of drag reduction we present the properties of Polyethy- 
lene oxide (N x [-CH2-CH2-O], known as PEO). The 
typical number N of monomers ranges between 10 4 and 
10 5 . In water solution in equilibrium this polymer is in 
a coiled state, with an end-to-end distance po of about 
1CP 7 m. When fully stretched the maximal end-to-end 
distance p m is about 5 x 10 -6 m. The typical mass con- 
centration used is between 10 to 1000 wppm. The viscos- 
ity of the water solution increases by a factor of 2 when 
PEO is added with concentration of 280 wppm. Note 
that even the maximal value value p m is much smaller 
than the typical Kolmogorov viscous scale, allowing us 
to model the turbulent velocity around such a polymer 
as a fluctuating homogeneous shear. 

Although the polymer includes many monomers and 
therefore many de grees of freedom, it was shown during 
years of research (IBeris and EdwardsL 11994 iBird et al, 
ll987tlde-Genned . ll979l : lFlorvlll953D that the most impor- 
tant degree of freedom is the end-to-end distance, allow- 
ing the simplest model of the polymer to be a dumbbell 
of two spheres of radius 5, and negligible mass, connected 
by a spring of equilibrium length po, characterized by 
a spring constant k. This model allows an easy under- 
standing of the coil-stretch transition under a turbulent 
shear flow. If such a dumbbell is stretched to length p p , 
the restoring force k(p p — po) is balanced by the Stokes 
force QiravQclpp/ dt (recall that the fluid density was taken 
as unity) . Thus the relaxation time r of the dumbbell is 
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r = Snauo/k. On the other hand, in the presence of a 
homogeneous shear S, the spring is stretched by a Stokes 
force ftiravoSpp. Balancing the stretching force with the 
restoring force (both proportional to p p for p p po), 
one finds that the coil-stretch transition is expected when 
St > 1. In fact, a free polymer will rotate under a homo- 
geneous shear, making this argument a bit more involved. 
Indeed, one needs a fluctuating shear s in order to stretch 
the polymer, leading finally to the condition 



velocity U a (r,t) as an additional stress tensor: 



> 1 . 



(3.1) 



To obtain a consistent hydrodynamic description of 
of polymer laden flows one needs to consider a field of 
polymers instead of a single chain. In a turbulent flow 
one can assume that the concentration of the polymers 
is well mixed, and approximately homogeneous. Each 
chain still has many degrees of freedom. However, as 
a consequence of the fact that the most important de- 
gree of freedom for a single chain is the end-to-end dis- 
tance, the effects of the ensemble of polymers enters the 
hydrodynamics in the form of a "conformation tensor" 
Rij (r, t) which stems from the ensemble average of the 
dyadic prod uct of the end-to-end dis t ance of the poly - 
mer chains (|Beris and Edwardsl 1 1994 iBird et all . 119871) . 
R a p = {p a pp). A successful model that had been em- 
ployed frequently in numerical simulations o f turbulent 
chan nel flows is the so-called FENE-P model ( Bird et al 
Il987t ). which takes into account the finite extensibility of 
the polymers: 



dR. 



a/3 



dt 



P(r,t) 



dU a 

(C/ 7 V 7 )i? Q /3 = ~q~~ RlP - 

-- [P(r,t)R a p - pl5 a(3 ] 

= {Pm~ Po)/(Pm~ Rjj) 



R 



(3.2) 
(3.3) 



The finite extensibility is reflected by the Peterlin func- 
tion P(r, t) which can be understood as a ensemble av- 
eraging correction to the potential energy k(p — po) 2 /2 
of the individual springs. For hydrodynamicists this 
equation should be evident: think about magneto- 
hydrodynamics, and the equations for the magnetic field 
n. Write down the equations for the diadic product 
n Q n^; the inertial terms will be precisely the first line 
of Eq. (|3.2p . Of course, the dynamo effect would then 
tend to increase the magnetic field, potentially without 
limit; the role of the second line in the equation is to 
guarantee that the finite polymer will not stretch with- 
out limit, and the Peterlin term guarantees that when 
the polymer stretches close to the maximum there will 
be rapid exponential decay back to equilibrium values of 
the trace of R. Since in most applications p m ^> po the 
Peterlin function can also be written approximately as 
P(r,t) w (1/(1 — 7-R 77 ) where 7 = p~ 2 . In its turn the 
conformation tensor appears in the equations for fluid 



dU a 
dt 



-(C/ 7 V 7 )C/a = -V Q p + v V 2 U a + V 7 T Q7 ,(3.4) 



T a p(r, t) 



P(r, t) 

5 Ra/3(r, t) - 8 a /3 

Po 



(3.5) 



Here v p is a viscosity parameter which is related to the 
concentration of the polymer, i.e. v p /vq ~ c p where c p 
is the volume fraction of the polymer. We note however 
that the tensor field can be rescaled to get rid of the pa- 
rameter 7 in the Peterlin function, R a p = ^Rap with 
the only consequence of rescaling the parameter v p ac- 
cordingly. Thus the actual value of the concentration is 
open to calibration against the experimental data. Also, 
in most numerical simulations, the term Pp~ 2 R is much 
larger than the unity tensors in (|3.2|) and (|3.5p . There- 
fore, in the theoretical development below we shall use 
the approximation 



T af 3 — v p PR a i3/Tp 



(3.6) 



We note that the conformation tensor always appears 
rescaled by p\. For notational simplicity we will absorb 
Pq into the definition of the conformation tensor, and 
keep the notation i? Q( g for the renormalized, dimension- 
less tensor. 

Considering first homogeneous and isotropic turbu- 
lence, The FENE-P model can be used to demonstrate 
the Lumley scale discussed above. In homogeneous tur- 
bulence one measures the moments of the velocity differ- 
ence across a length scale r', and in particular the second 
order structure function 

S 2 (r') = ({[U(r + r>, t) U(r, t)\ ■ ^} 2 ) , (3.7) 

r 

where the average (...) is performed over space and 
time. This quantity changes when the simulations 
are done with the Navier-Stokes equations on the one 
hand and with th e FEN E-P equations on the other, see 
( De-Angelis et all l2005j )). The polymers decrease the 
amount of energy at small scales. The reduction starts 
exactly at the theoretical estimate of the Lumley scale, 
shown with an arrow in Fig. El A more detailed analy- 
sis of the energy transfer ( De-Angelis et all . 120051 ) shows 
that energy flows from large to small scales and at the 
Lumley scale a significant amount of energy is transferred 
to the potential energy of the polymers by increasing i? 77 . 
This energy is eventually dissipated when the polymer re- 
lax their length back to e quilibrium. I n con tradistinction 
to the picture offered by Ide-Genned ( 1990f ). simulations 
indicate that the energy never goes back from the poly- 
mers to the flow; the only thing that the polymers can do 
is to increase the dissipation. This is a crucial statement 
underlying again the challenge of developing a consistent 
theory of drag reduction. In isotropic and homogeneous 
conditions the effect of polymer is just to lower the effec- 
tive Reynolds number which, of course, cannot explain 
drag reduction. 
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FIG. 6 Second order structure function S%(r) for homoge- 
neous and isotropic turbulence. Two cases are represented: 
with polymers (line with symbols) and without polymer 
(dashed line). The Lumley scale is indicated by an arrow. 
The numerical simulations with polymers is performed using 
a 128 3 resolution with periodic boundary conditions. The 
energy content for scales below the Lumley scale is reduced 
indicating a significant energy transfer from the velocity field 
to the polymer elastic energy. 



20 
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FIG. 7 Mean velocity profiles for the Newto nian and for 
the v iscoelastic simulations with Tie — 125 (Ben zi et all 
I2006T ). Solid line: Newtonian. Dashed line: Viscoelastic. 
The straight lines represent a log-law with the classical von 
Karman slope. Notice that in this simulation the modest 
Reynolds number results in an elastic layer in the region 
V + < 20. 



The FENE-P equations were also simulated on the 
computer in a channel or pipe geometry, reproduc- 
ing the phenomenon of drag reduction found in exper- 
imen t s (iBenzi et all 120061 : iDimitropoulos et all Il998l 
120051 iPtasinski et all l200lh . The most basic characteris- 
tic of the phenomenon is the increase of fluid throughput 
in the channel for the same pressure head, compared to 
the Newtonian flow. This phenomenon is demonstrated 
in Fig. [7J As one can see the simulation is limited 
compared to experiments; the Reynolds number is rel- 
atively low, and the MDR is not attained. Nevertheless 
the phenomenon is there. At any rate, once we are con- 
vinced that the FENE-P model exhibits drag reduction, 
it must also reproduce the universal properties of the 
phenomenon, and in particular the MDR. We will show 
that this is indeed the case, but that also all the cross- 
over non-universal phenomena can be understood using 
this model. If we were not interested in the cross-over 
phenomena we could use directly the large concentration 
limit of the FENE-P model, which is the limit P = 1 
where the model identifies with the harmonic Dumbbell- 
model known also as the Oldroid-B model. 



B. Derivation of the MDR 

In this section we review the theory that shows how 
the new log-law (jl.6p comes about when the polymers 
are added to the flow. In the next section we explain 
why this law is universal and estimate the parameters 
from first principles. 

As before, we proceed by taking the long time average 



of Eq. (|3.4[) , and integrating the resulting equation along 
the y coordinate. This produces an exact equation for the 
momentum balance: 

W + voS+^{PR xv )(y)=p'{L-y) . (3.8) 

T 

The interpretation of this equation is as before, but we 
have a new term which is the rate at which momentum 
is transferred to the polymers. Near the wall it is again 
permissible to neglect the term p'y on the RHS for y <C L. 
One should not be surprised with the form of Eq. (|3.8p . 
the new term could only be the one that is appearing 
there, since it must have ani-j/ symmetry, and it simply 
stems from the additional stress tensor appearing in Eq. 
p.4[) . The derivation of the energy balance equation is 
more involved, and had been described in great detail in 
(|Benzi et all 120061: iL'vov et all l2005afh The final form 
of the equation is 

K K 3 / 2 K(y) 
av a -+b +C41/ p {Ry )-m = W S, (3.9) 

yZ y y yZ 

where C4 is a dimensionless coefficient of the order of 
unity. This equation is in its final form, ready to be an- 
alyzed further. Equation (|3 . 8|) needs to be specialized to 
the vicinity of the MDR, which is only obtained when 
the concentration c p of the polymers is sufficiently high, 
when "Re is sufficiently high, but also when the D e borah 
number T>e is suffi ciently high. In (|Benzi et al. 1, 120061 
IL'vov et al.l . r2005al) it was shown that when these condi- 
tions are met, 

(PR xy ) = Cl (R yy )ST , (3.10) 
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where c\ is another dimensionless coefficient of the order 
of unity. Using this result in Eq. (|3.8[) we end up with 
the momentum balance equation 

W(y) + v S{y) + c lVp {R yy ) (y)S(y) = p'L , y^L . 

(3.11) 

The substitution (|3.10[) is so important for what follows 
that we dwell on it a bit further. When the conditions 
discussed above are all met, the polymers tend to be 
stretched and well aligned with the flow, such that the 
xx component of the conformation tensor must be much 
larger than the xy component, with the yy component 
being the smallest, tending to zero with T>e — > oo. Since 
De is the only dimensionless number that can relate the 
various components of (Rij), we expect that {R xx ) ~ 
Ve(R TM ) and (R xll ) ~ Ve{R yy ). Indeed, a calculation 
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( L'vov et al ., 2005a) shows that in the limit De — ► oo the 
conformation tensor attains the following universal form: 
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FIG. 8 Comparison of the DNS data Ref. (jSibila and Baronl 
120021 ) for the mean profiles of R xx and 10 x R yy , the compo- 
nents of the dimensionless conformation tensor, with analyt- 



forTte > 1 , (3.12) ica j prec iictions. In our notation = vapRL^Rij / (r v Rl 



where C is of the order of unity. We conclude that R xx 
is much larger than any other component of the confor- 
mation tensor, but it plays no direct role in the phe- 
nomenon of drag reduction. Rather, R yy which is very 
much smaller, is the most important component from 
our point of view. We can thus rewrite the two balance 
equations derived here in the suggestive form: 



v(y)S + W = p'L , 

K K 3 / 2 

av{y)^+b = WS 

y y 



(3.13) 
(3.14) 



where a — ac^jc\ and the "effective viscosity" v(y) is 
being identified as 



v(y) = v + civ p (R yy ) 



(3.15) 



Exactly like in the Newtonian theory one needs to add 
a phenomenological relation between W(y) and K(y) 
which holds in the elastic layer, 



W(y) = c 2 K(y) 



(3.16) 



with c v an unknown coefficient of the order of unity. 

At this point one asserts that the 'dressed' viscous term 
dominates the inertial term on the LHS of the balance 
equations (|3.13p and (|3.14|) . From the first of these we 
estimate v(y) ~ p'L/S(y). Plugged into the second of 
these equations this leads, together with Eq. (|2.6[) to 
S(y) =Const./y, where Const, is a combination of the 
unknown coefficients appearing in these equations, and 
therefore is itself an unknown coefficient of the order of 
unity. The important thing is that the theory predicts a 
new log-law for V(y), the slope of which we will shown to 
be universal in the next Subsection. Before we do so it is 
important to realize that if S(y) =Const./y, our analysis 



Black squares: DNS data for the streamwise diagonal compo- 
nent R xx , that according to our theory has to decrease as 1 /y 
with the distance to the wall. Black solid line: the function 
l/y + ■ Red empty circles: DNS data for the wall-normal com- 
ponent WRyy, for which we predicted a linear increase with 
y + in the log-law turbulent region. Red dashed line - linear 
dependence, oc y + . 



indicates that the effective viscosity u(y) must grow out- 
side the viscous layer linearly in y, v{y) ~ y. If we make 
the self-consistent assumption that i/q is negligible in the 
log-law region compared to v p (R yy ), then this prediction 
appears in contradiction with what is known about the 
stretching of polymers in channel geometry, where it had 
been measured that the extent of stretching decreases as 
a function of the distance from the wall. The apparent 
contradiction evaporates when we recall that the amount 
of stretching is dominated by R xx which is indeed de- 
creasing. To see this note that Eq. (|3.12[) predicts that 
R xx ~ 2T>e 2 R yy (x 1/y. At the same time R yy increases 
linearly in y. Both predictions agree with what is ob- 
served in simulations, see Fig. [5] 

The simplicity of the resulting theory, and the cor- 
relation between a linear viscosity profile and the phe- 
nomenon of drag reduction raises the natural question, 
is it enough to have a viscosity that rises linearly as 
the function of the distance from the wall to cause 
drag reduction? To answer this question one can sim- 
ulate the Navier-Stokes equations with proper viscos- 
ity profiles (discussed below) and show that the re- 
sults are in semi-quantitative agreement with the corre- 
sponding full FENE-P D NS. Such simulations were done 
(|De-Angelis et all \2004 in a domain 2ttL x 2L x 1.2-kL, 
with periodic boundary conditions in the streamwise and 
spanwise directions, and with no slip conditions on the 
walls that were separated by 2L in the wall-normal di- 
rection. An imposed mass flux and the same Newtonian 
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FIG. 9 The Newtonian viscosity profile and four examples 
of close to linear viscosity profiles employed in the numerical 

simulations. Solid line: run N, : run R, ■ • •: run S, — ■ — : 

T, : run U. 




FIG. 10 The reduced mean velocity as a function of the re- 
duced distance from the wall. The line types correspond to 
those used in Fig. [9] 



initial conditions were used. The Reynolds number Re 
(computed with the centerline velocity) was 6000 in all 
the runs. The y dependence of the scalar effective viscos- 
ity was close to being piece- wise linear along the channel 
height, namely v = v$ for y < j/i, a linear portion with 
a prescribed slope for y\ < y < j/2, and again a constant 
value for y2 < y < L. For numerical stability this profile 
was smoothed out as shown in Fig. [9] Included in the fig- 
ure is the flat viscosity profile of the standard Newtonian 
flow. In Fig. [10] we show the resulting profiles of V^~{y) 
vs. y + . The line types are chosen to correspond to those 
used in Fig. [5] The decrease of the drag with the increase 
of the slope of the viscosity profiles is obvious. Since the 
slopes of the viscosity profiles are smaller than needed to 
achieve the MDR asymptote for the corresponding 7?e, 
the drag reduction occurs only in the near-wall region 
and the Newtonian plugs are clearly visible. 

The conclusion of these simulations appears to be that 
one can increase the slope of the linear viscosity profile 
and gain further drag reduction. The natural question 
that comes to mind is whether this can be done without 
limit such as to reduce the drag to zero. Of course this 
would not be possible, and here lies the clue for under- 
standing the universality of the MDR. 



C. The universality of the MDR: theory 

The crucial new insight that will explain the universal- 
ity of the MDR and furnish the basis for its calculation 
is that the MDR is a marginal flow state of wall- 
bounded turbulence: attempting to increase S(y) be- 
yond the MDR results in the collapse of the turbulent so- 
lutions in favor of a stable laminar solution with W — 
( Benzi et al l l200fiaj> . As such, the MDR is universal by 
definition, and the only question is whether a polymer 
(or other additive) can supply the particular effective vis- 



cosity v{y) that drives Eqs. (|3.13[) and (|3.14| to attain 
the marginal solution that maximizes the velocity pro- 
file. We expect that the same marginal state will exist in 
numerical solutions of the Navier-Stokes equations fur- 
nished with a y-dependent viscosity v(y). There will be 
no turbulent solutions with velocity profiles higher than 
the MDR. 

To see this explicitly, we first rewrite the balance equa- 
tions in wall units. For constant viscosity [i.e. v(y) = vq], 
Eqs. (|3.13[) - (|2.6jl form a closed set of equations for 
S + = Sv /(p'L) and W + = W/p'L in terms of two di- 
mensionless constant S + = a^J K/W (the thickness of 
the viscous boundary layer) and k k = b/c^ (the Von 
Karman constant). Newtonian experiments and simula- 
tions agree well with a fit using S + ~ 6 and k k ~ 0.436 
(see the black continuous line in Fig. [2] which shows the 
mean velocity profile using these very constants). Once 
the effective viscosity v{y) is no longer constant we ex- 
pect c N to change (c N — * c v ) and consequently the two 
dimensionless constants will change as well. We will de- 
note the new constants as A and k c respectively. Clearly 
one must require that for v(y)/vo — > 1, A — > S + and 
The ba lance equations are now written as 
(|Benzi et al.l . l2005ah : 

v + (y + )S+(y+) + W+(y+) = 1 , (3.17) 

^+)^ + ^ = S + , (3.18) 

y + K cV + 

where v + {y + ) = v{y + )/v . Substituting now S + from 
Eq. (|3.17ll into Eq. (|3.18p leads to a quadratic equation 
for VW + . This equation has as a zero solution for W + 
(laminar solution) as long as v + (y + )A/y + = 1. Turbu- 
lent solutions are possible only when v + (y + )/S./y + < 1. 
Thus at the edge of existence of turbulent solutions we 
find u + ex y + for y + 3> 1. This is not surprising, since it 
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was observed above that the MDR solution is consistent 
with an effective viscosity which is asymptotically linear 
in y + . It is therefore sufficient to seek the edge solu- 
tion of the velocity profile with respect to linear viscosity 
profiles, and we rewrite Eqs. (|3.17[) and (|3.18[) with an 
effective viscosity that depends linearly on y + outside the 
boundary layer of thickness S + : 



[1 + a(y+ - 5+)}S+ + W+ = 1 , 



5 H 



(3.19) 
(3.20) 



We now endow A with an explicit dependence on the 
slope of the effective viscosity v + {y), A = A(a). Since 
drag reduction must involve a decrease in W, we expect 
the ratio a 2 K/W to depend on a, with the constraint 
that A (a) — > S + when a — > 0. Although A, <5+ and 
a are all dimensionless quantities, physically A and S + 
represent (viscous) length scales (for the linear viscos- 
ity profile and for the Newtonian case respectively) while 
or 1 is the scale associated to the slope of the linear vis- 
cosity profile. It follows that a<5 + is dimensionless even in 
the original physical units. It is thus natural to present 
A (a) in terms of a dimensionless scaling function f(x), 



A(a) = 6+f{a5+) 



(3.21) 



Obviously, /(0) = 1. In the Appendix we show that the 
balance equations (]3 . 19[) and (|3.20|) (with the prescribed 
form of the effective viscosity profile) have an non-trivial 
symmetry that leaves them invariant under rescaling of 
the wall units. This symmetry dictates the function A(a) 
in the form 



A(a) = 



6+ 



1 - aS+ 



(3.22) 



Armed with this knowledge we can now find the maximal 
possible velocity far away from the wall, y + 3> S + . There 
the balance equations simplify to 



ay+S + + W + = 1 , 
aA 2 {a) + VW+/n c = y+S^ 



(3.23) 
(3.24) 



T hese equations have the j/ + -independent solution for 
VW+ and y+S+: 



y + S + = aA 2 (a) + VW+/k c . (3.25) 



By using equation (|3.25ll , (see Fig. QT]), we obtain 
that the edge solution (W + — > 0) corresponds to the 
supremum of y + S + , which happens precisely when a — 
1/A(a). Using Eq. (|3.22[) we find the solution a — a m — 
1/2(5+. Then y+S+ = A(a m ), giving k' 1 = 25+. Using 
the estimate 5+ « 6 we get the final prediction for the 
MDR. Using Eq. CEI|) with k" 1 = 12, we get 
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FIG. 11 The solution for WVW+ (dashed line) and y+S + 
(solid line) in the asymptotic region y + S + , as a function 
of a. The vertical solid line a = l/25 + = 1/12 which is the 
edge of turbulent solutions; Since VW+ changes sign here, 
to the right of this line there are only laminar states. The 
horizontal solid line indicates the highest attainable value of 
the slope of the MDR logarithmic law 1/k v = 12. 



V + {y + ) w 121ny+ - 17.: 



(3.26) 



This result is in close agreement with the empirical law 
(|1.5j) proposed by Virk. The value of the intercept on the 
RHS of Eq. (pOtJj) follows from Eq. (JTBJ) which is based 
on matching the vi scous solution to the MDR log-law in 
( L'vov et all I2004T ). We now also have the deep justifi- 
cation for this matching: the MDR is basically a laminar 
solution that can match smoothly with the viscous sub- 
layer, with continuous derivative. This is not possible for 
the von Karman log-law which represents fully turbulent 
solutions. Note that the numbers appearing in Virk's law 
correspond to 6 + = 5.85, which is well within the error- 
bar on the value of this Newtonian parameter. Note that 
we can easily predict where the asymptotic law turns into 
the viscous layer upon the approach to the wall. We can 
consider an infinitesimal W + and solve Eqs. (|3.17p and 
p,18p for S + and the viscosity profile. The result, as be- 
fore, is v + {y) = A{a m )y + . Since the effective viscosity 
cannot fall bellow the Newtonian limit v + = 1 we see 
that the MDR cannot go below y+ = A(a m ) = 28+ . We 
thus expect an extension of the viscous layer by a factor 
of 2, in very good agreement with the experimental data. 

Note that the result W + = should not be interpreted 
as W — 0. The difference between the two objects is the 
factor of p'L, W — p'LW + . Since the MDR is reached 
asymptotically as 7?e — > oo, there is enough turbulence 
at this state to stretch the polymers to supply the needed 
effective viscosity. Indeed our discussion is in close corre- 
spon dence with the experimental remark by Virk ( Virkl . 
Il975l ) that close to the MDR asymptote the flow appears 
laminar. 

In summary, the main message of this section is that 
the added polymers endow the fluid with an effective vis- 
cosity v(y) instead of v§. There exists a profile of v(y) 
that results in a maximal possible velocity profile at the 
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edge of existence of turbulent solutions. That profile is 
the prediction for the MDR. In particular we offer a pre- 
diction for simulations: direct numerical simulations of 
the Navier-Stokes equations in wall bo unded geometries , 
endow ed with a linear viscosity profile (|De-Angelis et all . 
2004). will not be able to support turbulent solutions 
when the slope of the viscosity profile exceeds the criti- 
cal value that is in correspondence with the slope of the 
MDR. 



IV. THE "ADDITIVE EQUIVALENCE": THE MDR OF 
ROD-LIKE POLYMERS 

In this section we address the experimental finding that 
rigid rod-like polymers app ear to exhibit the same MDR 
([T6]) as flexible polymers (|Virk et all Il997l ). Since the 
bare equations of motion of rod-like polymers differ quite 
significantly from those of flexible polymers, one needs to 
examine the issue carefully to understand this universal- 
ity, which was termed by Virk "Additive Equivalence" . 
We will see that the point is that in spite of the different 
basic equations, when the conditions allow attainment of 
the MDR, the balance equations for momentum and en- 
ergy are iden t ical in form to those of the flexible polymers 
(|Benzi et al.[ l2005bl : IChTng et all . [20061 ). The differences 
between the two types of polymers arise when we consider 
how the MDR is approached, and cross-overs back to the 
Newtonian plugs, all issues that are taken up below. 



where v p is the polymeric contribution to the viscosity 
at vanishingly small and time-independent shear; v p in- 
creases linearly with the polymer concentration, making 
it an appropriate measure for the polymer's concentra- 
tion. The other quantities in Eq. (|4. 2[) are the velocity 
gradient tensor 



S ab = dU a /dx b , 



(4.3) 



and n = n(r, t) is a unit (n-n = 1) director field that de- 
scribe the polymer's orientation. Notice, that for flexible 
polymers the equation (|3.6p for T a b is completely different 
from (|4~2]) . The difference between Eqs. ((472]) and ([376]) 
for the additional stress tensor in the cases of rod-like and 
flexible polymers reflects their very different microscopic 
dynamics. For the flexible polymers the main source of 
interaction with the turbulent fluctuations is the stretch- 
ing of the polymers by the fluctuating shear s. This is 
how energy is taken from the turbulent field, introducing 
an additional channel of dissipation without necessarily 
increasing the local gradient. In the rod- like polymer case 
the dissipation is only taken as the skin friction along the 
rod-like polymers. Having in mind all these differences 
it becomes even more astonishing that the macroscopic 
equations for the mechanical momentum and kinetic en- 
ergy balances are isomorphic for the rod-like and flexible 
polymers, as is demonstrated below. 



B. The balance equations and the MDR 



A. Hydrodynamics with rod-like polymers 

The equation for the incompressible velocity field 
U(r,t) in the presence of rod- like polymers has a form 
isomorphic to Eq. (|3.4[) 

+ VU = v &U- Vp + V-a , (4.1) 

but with another cr =>■ a a b playing the role of an extra 
stress tensor caused by the polymers. 

The calculation of t he tensor a for r i gid ro ds is of- 
fered in the literature ( Doi and Edward si Il988l ). subject 
to the realistic assumptions that the rod-like polymers 
are mass-less and having no inertia. In other words, the 
rod-like polymers are assumed to be at all times in local 
rotational equilibrium with the velocity field. Thus the 
stress tensor does not have a contribution from the rota- 
tional fluctuations against the fluid, but rather only from 
the velocity variations along the rod-like object. Such 
variations lead to "skin friction" , and this is th e only ex- 
tra d i ssipative effect that is t aken i nto account dBrennerl . 
119741: iHinch and Leal . Il975l 119761 : iManhartl . |2003j). The 
result of these considerations is the following expression 
for the additional stress tensor: 

o a b = 6u p n a nb {niUjSij) , rod-like polymers , (4.2) 



Using Eq. (|4.2p and Reynolds decomposition (|2.1[) we 
compute 

(<Txy) = QypiUxyTlijSij) = 6v p [S(H? ) + (TZxyTZijSij)} . 

(4.4) 

Now we make use of the expected solution for the con- 
formation tensor in the case of large mean shear. In 
such flows we expect a strong alignment of the rod-like 
polymers along the streamwise direction x. The direc- 
tor components n y and n z are then much smaller than 
n x w 1. For large shear we can expand n x according to 



(4.5) 



This expansion allows us to express all products 
TZ-abTZ-cd = n a nbn c rid in terms that are linear in 7t, up 
to third order terms in n y ~ n z . In particular 



ni 



1 - 2(Ky 



V 2 7? 
,x -xy ,x, yy 



(4.6) 



With (TQ|) the first term in the RHS of Eq. (|4~4"]l can be 

estimated as 



6V p S(TZl ) = Clfp(Ryy)S, 



ci ~ 6 



(4.7) 



The estimate of the second term on the RHS of Eq. 
needs some f urther calculations which can be found in 
(jBenzi et all l2005bl : IChing et all 120061 ) . with the result 
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that it is of the same order as the first one. Finally we 
can present therefore the momentum balance equation in 
the form 



v S + c 1 v p {R yy )S + W = p'L 



(4.8) 



Another way of writing this result is in the form of an 
effective viscosity, 



v(y)S + W=p'L , 



(4.9) 



where the effect of the rod-like polymers is included by 
the effective viscosity v{y): 



v{y) = v + civ p (R yy ) 



(4.10) 



We see that despite the very different microscopic stress 
tensors, the final momentum balance equation is the 
same for the flexible and the rod-like poly mers. Addi- 
tional calculations which can be found in (|Benzi et all 
2005b) show that the energy balance equation also at- 
tains the precisely the same form as Eq. (|3.14[) . Evi- 
dently this immediately translates, via the theory of the 
previous sections, to the same MDR by the same mech- 
anism, and therefore the "additive equivalence" . 



V. NON-UNIVERSAL ASPECTS OF DRAG REDUCTION: 
FLEXIBLE POLYMERS 

In this section we return to the cross-over phenom- 
ena described in Subsect. II.B.2I and Su bsect. II.B.3I and 
provide the theory ( Benzi et all l2004bh for their under- 
standing. First we refer to the experimental data in Fig. 



1 



A. The efficiency of drag reduction for flexible polymers 

When there exist cross-overs back from the MDR to 
the Newtonian plugs the mean velocity pr ofile in the fl ex- 
ible polymer case consists of three regions (IVirkl . ll975h : a 
viscous sublayer, a logarithmic elastic sublayer (region 3 
in the Fig. [T2|) with the slope greater then the Newtonian 
one, Eq. II. 5p . and a Newtonian plug (region 4). In the 
last region the velocity follows a log-law with the New- 
tonian slope, but with some velocity increment AV + : 



V + (y+) = K- 1 \ny+ + B K +AV J 



(5.1) 



Note that we have simplified the diagram for the sake of 
this discussion: the three profiles Eqs. (|1.2|) , (|1.3[) , and 
(II. 5|) intercept at one point y + = y+ ~ n^ 1 ~ 11.7 ~ 
2<5 + . In reality the Newtonian log-law does not connect 
sharply with the viscous solution V + = y + . but rather 
through a cross-over region of the order of 5 + . 

The increment AV + which determines the amount of 
drag reduction is in turn determined by the cross-over 
from the MDR to the Newtonian plug (see Fig. [T2"|) . We 



)>: y + y log y + 

FIG. 12 Schematic mean velocity profiles. Region 1: y + < 
y J, - viscous sublayer. Region 2: y + > y+ - logarithmic 
layer for the turbulent Newtonian flow. Region 3: y^ < y + < 
y+ = (l + Q)^ 
flow. Region 4: 
flow. 



MDR asymptotic profile in the viscoelastic 
Newtonian plug in the viscoelastic 



y + > vt 



refer to this cross-over point as y+ . To measure the qual - 
ity of drag reduction one introduces ( Benzi et all l2004bh 
a dimensionless drag reduction parameter 



Q 



(5.2) 



The velocity increment AV + is related to this parameter 
as follows 

AV+ = (k- 1 - k- 1 ) In (y+/y+) = /31n(l + Q) . (5.3) 

Here f3 = k" 1 — k" 1 ~ 9.4. The Newtonian flow is then 
a limiting case of the viscoelastic flow corresponding to 
Q = 0. 

The cross-over point z/~J~ is non-universal, depending 
on Tie, on the number of polymers per unit volume c p , 
the chemical nature of the polymer, etc. According to 
the theory of the last Section, the total viscosity of the 
fluid v to t(y + ) = v + v p {y + ) [where v p {y + ) is the poly- 
meric contribution to the viscosity which is proportional 
to (Ryy)} is linear in y + in the MDR region: 



^tot(?/ + ) = vvy + /yl 



vZ <y < vZ 



(5.4) 



When the concentration of polymers is small and Tie is 
large enough, the cross-over to the Newtonian plug at j/"J" 
occurs when the polymer stretching can no longer provide 
the necessary increase of the total fluid viscosity. In other 
words, in that limit the cross-over is due to the finite 
extensibility of the polymer molecules. Obviously, the 
polymeric viscosity can not be greater than v v max which 
is the viscosity of the fully stretched polymers. Thus 
the total viscosity is limited by vq + P pm ax- Equating 
v p max and vtotiy*) gives us the cross-over position 



(5.5) 
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It follows from Eq. (|5-2[) that the drag reduction param- 
eter is determined very simply by 



Q 



y p max 



/vq , c p small, T^e large 



(5.6) 



At this point we need to relate the maximum poly- 
meric viscosity ^ pma x to the polymer properties. To this 
aim we estimate the energy dissipation due to a single, 
fully stretched, polymer molecule. In a reference frame 
co-moving with the polymer's center of mass the fluid ve- 
locity can be estimated as u ~ rVu (the polymer's center 
of the mass moves with the fluid velocity due to negli- 
gible inertia of the molecule) . The friction force exerted 
on the «-th monomer is estimated using Stokes law with 
Sui being the velocity difference across a monomer, 



p u a Sui = PqVq a r,V« , 



(5.7) 



where a is an effective hydrodynamic radius of one 
monomer (depending on the chemical composition), and 
r,; is the distance of the i-th monomer from the center 
of the mass. In a fully stretched state ~ ai (the 
monomers are aligned along a line). The energy dissi- 
pation rate (per unit volume) is equal to the work per- 
formed by the external flow 



dE 
~dt 



N p 

i=l 
PQ^pmax 



/9 iWcpA^(Vw) 2 



We thus can estimate v. 



p max • 
Vpma.x = Voa 3 C p Np 



(5.8) 



(5.9) 



Finally, the drag reduction parameter Q is given by 

Q = a 3 c p Np c p small, 7?e large . (5.10) 
This is the central theoretical result of iBenzi et al.l 



(2004b), relating the concentration c p and degree of poly- 
merization N p to the increment in mean velocity AV" 1 
via Eq. (15. 



B. Drag reduction when polymers are degraded 

The main experimental results which are of interest 
to us are summarized in Fig. [3J Note that in this ex- 
periment the flow geometry is rather complicated: with 
counter-rotating disks the linear velocity depends on the 
radius, and the local Reynolds number is a function of 
the radius. The drag reduction occurs however in a rela- 
tively small near-wall region, where the flow can be con- 
sidered as_a L jlow_neax_the flat plate. Thus, one con- 
siders ( Benzi et all l2004b[ ) an equivalent channel flow - 
with the same Tie and a half width L of the order of 
height/radius of the cylinder. In this plane geometry the 
torques in (| 1 . Tj) should be replaced by the pressure gra- 
dients p' : 



%DR: 



Pi 



x 100 



(5.11) 




4 6 
DR parameter a 

FIG. 13 Drag reduction efficiency %DR as function of the 
drag reducing parameter Q for different centerline Reynolds 
numbers Re: 1.2 x 10 s , 1.2 x 10 6 , and 1.2 x 10 7 (from top to 
bottom). 



In order to relate %DR with the drag reduction param- 
eter Q, one re- writes Eq. (15. ip in natural units 



(5.12) 

To find the degree of drag reduction one computes p' v and 
p' n keeping the centerline velocity V(L) constant. Defin- 
ing the centerline Reynolds number as Re, we rewrite 

Re = V ( L ) L = -fo, \ K -i ln7 fe + B K + AF+] . (5.13) 

This equation implicitly determines the pressure gradi- 
ent and therefore the %DR as a function of Q and Re. 
The set of Eqs. (15. 3[) and (|5.13[) is readily solved nu- 
merically, and the solution for three different values of 
Re is shown in Fig. 1131 The middle curve corresponds 
to Re= 1.2 x 10 6 , which coinc ides with the experimental 



conditions (Choi et al ., 2002). One sees, however, that 
the dependence of %DR on Re is rather weak. One im- 
portant consequence of the solutions shown in Fig. 1131 is 
that for small Q (actually for Q < 0.5 or %DR < 20), 
%DR is ap proximately a lin ear function of Q. The ex- 
periments ( Choi et al. , [2002) lie entirely within this lin- 
ear regime, in which we can linearize Eq. (|5 . 13[) in Ay + , 
solve once for p' v and once for p' K (using AV + = 0). 
Computing Eq. I|5.11[) we find an approximate solution 
for the %DR: 



%DR 



2/3Q 



■ln(7& N )+B K 



x 100 



(5.14) 



Here 7Ze N is the friction Reynolds number for the New- 
tonian flow, i.e. the solution of Eq. (|5 . 13|) for AV + = 0. 

It is interesting to note, that while the %DR depends 
on the Reynolds number, the ratio of different %DR's 
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does not [to 0(Q)\. 

%DR ( 



Q 



(i) 



%DR 



(2) 



i/p max 



(5.15) 



This result, together with Eq. ()5.9j) , rationalize s com - 
pletely the experimental finding of ( Choi et"aT1 . |2002| ) 
summarized in Fig. [31 During the DNA degradation, 
the concentration of polymers increases by a factor of 2, 
while the number of monomers N p decreases by the same 
factor. This means that %DR should decrease by a factor 
of 4, as is indeed the case. 

The experimental results pertain to high Tie and small 
c p , where we can assert that the cross-over results from 
exhausting the stretching of the polymers such that the 
maximal available viscosity is achieved. In the linear 
regime that pertains to this experiment the degradation 
has a maximal effect on the quality of drag reduction Q, 
leading to the precise factor of 4 in the results shown in 
Fig. [3J Larger values of the concentration of DNA will 
exceed the linear regime as is predicted by Fig. [T31 then 
the degradation is expected to have a smaller influence 
on the drag reduction efficacy. It is worthwhile to test the 
predictions of this theory also in the nonlinear regime. 



C. Other mechanism for cross-over 

Having any reasonable model of polymer-laden flows at 
our disposal we can address other possibilities for the sat- 
uration of drag reduction. We expect a cross-over from 
the MDR asymptote back to the Newtonian plug when 
the basic assumptions on the relative importance of the 
various terms in the balance equations lose their valid- 
ity, i.e. when (i) turbulent momentum flux W becomes 
comparable with the total momentum flux p' L, or when 
(ii) turbulent energy flux bK 3 / 2 /y becomes of the same 
order a s turbulent energy production WS. In fact it was 
shown (|Benzi et al.l . l2006t) using the FENE-P model that 
both these conditions give the same cross-over point 




(5.16) 



Note that f(y) = r/(P(y)) is the effective non-linear 
polymer relaxation time. Therefore condition (|5.16p can 
be also rewritten as 



(5.17) 



In writing this equation we use the fact that the cross- 
over point belongs also to the edge of the Newtonian 
plug where S(y) ~ \Jp'Ljy. The LHS of this equation is 
simply the local Deborah number (product of local mean 
shear and local effective polymer relaxation time). Thus, 
the cross-over to the Newtonian plug ocurrs at the point, 
where the local Deborah number decreases to ~ 1. Wc 
expect that this result is correct for any model of elastic 
polymers, not only for the FENE-P model considered 
here. 



To understand how the cross-over point y v depends 
on the polymer concentration and other parameters, one 
needs to estim ate mean value of the Peterlin function 
(P). Following iBenzi et al. I (120061) wc estimate the value 
of (P) as 

( p ) = ^-77^ > (5-18) 



Pi, 



1 - j(R) ' 

Rzz) ~ (Rxx 



'iRyy) and 



where (R) = (R xx -r ^ yy 

7 ~ (for simplicity we disregard po). We know 

from before that 



(5.19) 



(R xx ) ~ (Sf) 2 (R yy ) 
and therefore at the cross-over point (|5.17p 

{Rxx) - (Ryy) ■ (R) - (Ryy) ■ 

The dependence of (R yy ) on y in the MDR region follows 
from Eqs. (|3.11[) , since 




Then at the cross-over point y = y^ 



1 - IVvVWL/vp ' 
Substituting this estimate into (|5. 16[> gives the final result 



1 + ■yp'Lr/ v p 



(5.20) 



Here C is constant of the order of unity. Finally, intro- 
ducing dimensionless concentration of polymers 

. _ Vp 
Cp = , 

7^o 



(5.21) 



one can write denominator in Eq. (|5 . 20[) as 
jp' L T 



1 



where 



= 1 



Ve 



1 p'Lr Ve 
= 1 + — 

Zp Vq c v 



p'Lr 

V<3 



(5.22) 



is the (global) Deborah number. Then for the dimension- 
less cross-over point y"J~ = y w \fp'LfvQ one obtains 

CVe 



1 + De/c p 



(5.23) 



This prediction can be put to direct test when c p is 
very la rge, or equivalently in the Oldroyd B model where 
P = 1 (|Benzi et al.l . l2004al) . Indeed, in numerical simula- 
tions when t he Weissenberg number was changed system- 
atically, cf. (|Yu et al.l . [200l . one observes the cross-over 
to depend on TJe in a manner consistent with Eq. (|5.23[) . 
The other limit when c p is very small is in agreement with 
the linear dependence on c p predicted in Eq. (| 5 . 1 [) . 

We can thus reach conclusions about the saturation of 
drag reduction in various limits of the experimental con- 
ditions, in agreement with experiments and simulations. 
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FIG. 14 The mean velocity profiles for flexible polymer ad- 
ditives with v = 1,5,10,20,50,100 and 500 from below to 
above. Note that the profile follows the MDR until it crosses 
over back to the Newtonian plug. 



FIG. 15 The mean velocity profile for rod-like polymer ad- 
ditives with u = l,5, 10, 20, 50, 100, 500, 1000, 5000 and 10000 
from below to above. Note the typical behavior expected for 
rod-like polyemrs, i.e. that the profile diverges from the von 
Karman log-law, reaching the MDR only asymptotically. 



VI. CROSS-OVER PHENOMENA WITH ROD-LIKE 
POLYMERS 

A. The attainment of the MDR as a function of 
concentration 

At this point we return to the different way that the 
MDR is approached by flexible and rod-like polymers 
when the concentration of the polymer is increased. As 
discussed above, in the case of flexible polymers the MDR 
is followed until the cross-over point y+ that was dis- 
cussed already in great detail. The rod-like polymers 
attain the MDR only asymptotically, and for intermedi- 
ate values of the concentration the mean velocity profile 
increase gradually from the von Karman log-law to the 
MDR log-law. 

This difference can be fully understood in the context 
of the present theory. The det ailed calcula t ion a ddress- 
ing this issue was presented in ( Ching et all [2006) . The 
results of the calculation, presented as the mean velocity 
profiles for increasing concentration of the two types of 
polymers are shown in Figs. [TH and [15j The reader 
should note the difference between these profiles as a 
function of the polymer concentration . While the flex- 
ible polymer case exhibits the feature ( Virk et all . Il997t 
?) that the velocity profile adheres to the MDR until a 
cross-over to the Newtonian plug is realized, the rod-like 
case presents a "fan" of profiles which only asymptot- 
ically reach the MDR. We also notice that the flexible 
polymer matches the MDR faithfully for relatively low 
values of v = v p /vo, whereas the rod-like case attains the 
MDR only for much higher values of v. Thi s result is in 
agreement with the experi mental finding in ( Bonn et all 
2005|; IWagner et all 120031 ) that the flexible polymer is a 



better drag reducer than the rod-like a nalogu e. 

The calculation of (|Ching et all 120061) allows a 
parameter-free estimate of the cross-over points from 
the MDR to the "newtonian plug" in the case of flexible 
polymers. The resulting estimate reads 
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T).l . 



(6.1) 



These estimates agrees well with the numerical results in 
Fig. [14] No such simple calculation is available for the 
case of the rod-like polymers since there is no clear point 
of departure for small v. 

We should note that the higher efficacy of flexi- 
ble polymers cannot be easily related to their elonga- 
tional vi scosity as m e asured in laminar flows. I n some 
studies (iBonn et all 120051: Iden Toonder et all Il995t 
IWagner et all 20031 ) it was proposed that there is a cor- 
relation between the elongational viscosity measured in 
laminar flows and the drag reduction measured in turbu- 
lent flows. We find here that flexible polymers do better 
in turbulent flows due to their contribution to the ef- 
fective shear viscosity, and their improved capability in 
drag reduction stems simply from their ability to stretch, 
something that rod-like polymers cannot do. 



B. Cross-over phenomena as a function of the Reynolds 
number 

Finally, we address the drag enhancement by rod-like 
polymers when the values of Tie are t oo sm all, see Fig. [5] 
The strategy of lAmarouchene et al.l ( 2007| ) is to develop 
an approximate formula for the effective viscosity in the 
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case of rod-like polymers that interpolates properly be- 
tween low and high values of Tie. For that purpose we 
remind the reader of the high Tie form of v c s, derive the 
form for low value of Tie, and then offer an interpolation 
formula. 

Even for intermediate values of Tie one cannot neglect 
y in the production term p'(L-y) in the momentum bal- 
ance equation. Keeping this term the momentum balance 
equation becomes, in wall units, 



(1 + u+R w )S + + W + = l- a —, 



where 2/+ = v p /vq. 

The energy balance equation remains as before, 



(6.2) 



K+ (K + \ 3 / 2 
W+S + ~a\l + v+R yy )—+b [ —^- . (6.3) 



(y + ) 2 



As was explained above, equations (|6.2p and (|6.3p imply 
that the polymers change the properties of the flows by 
replacing the viscosity by 



^cff = 1 + v+R 



p -"-TO 



(6.4) 



In the fully developed turbulent flow with rod-like 
polymers, when Tie is very large, it was shown in 
( Benzi et all l2005bh that R yy depends on K + and S + : 



Ryy — 



K- 



{y+S+f 



(6.5) 



It was argued in ( Benzi et a l.. 2005b) that for large Re, 
K + grows linearly with y + and thus the viscosity profile 
is linear. 

Next consider low Tie flows. According to Eq. (|6.4p . 
the value of v c g depends on 1/+ and R yy . The value 
of 1/+ is determined by the polymer properties such as 
the number of monomers, their concentration etc., and 
thus 1/+ should be considered as an external parame- 
ter in the equation. The value of R yy , on the other 
hand, depends on the properties of the flow. In the 
case of laminar flow with a constant shear rate, i.e., 
K + = W + = and S + =constant, it was shown the- 
oretically in ( Doi and Edward's! Il988l ) that: 



2 i/3 



R yy - 0,2/3 



(6.6) 



Thus, the effective viscosity is reduced if S is increased, 
and therefore the rod-like polymers solution is a shear- 
thinning liquid. Naturally, The value of De changes with 
Tie. To clarify this dependence we consider the momen- 
tum equations Eq. (|3.13p at y = in the Newtonian case. 



and so Tie depends on p'L only. According to (jl.l|) . Tie T 
grows as \Jp' L and therefore 

^0 



Ve = 



7L 2 



Tie z 



by Eq. (|6.7p . Putting into Eq. (|6.6p . we have 



^eff 



1 



A 



P Tie 4 / 3 



(6.8) 



(6.9) 



where A = vq/^L 2 is a constant. 

In the case of intermediate Tie, we need an interpola- 
tion between Eqs. (|6.5p and (|6.6p . To do this we note 
that when y + is small, the solution of Eqs. (|2.6p . (|6.2p 
and (|6.3p result in W + = K + = in the viscous sub- 
layer. This implies that the flow cannot be highly tur- 
bulent in the viscous sub-layer. Thus, it is reasonable 
to employ Eq. (|6.6p as long as y + is small. On the other 
hand, as the upper bound of y + is Tie, when y + is large, it 
automatically implies that Tie. is large. The laminar con- 
tribution is therefore negligible as it varies inversely with 
Tie. The effective viscosity due to the polymer is domi- 
nated by the turbulent estimate, Eq. (|6 . 5[) . To connect 
these two regions we simply use the pseudo-sum: 



z/ cff = 1 + 



X 



K 4 



.9 + K 



p \Tie 4 / 3 {y+S+) 2 
K+ 



(y+S+f 



(6.10) 



where g = 1 + v+ X/Tie 4 ^ 3 . One can see that the limits 
for both high and low Tieare satisfied. 

The form for (|6.10p for i/ e g was incorporated into the 
balance equatio ns which were analyzed and solved self- 
consistently in ( Amarouchene et all l2007f ) one of the 
most interesting results, which can be directly compared 
to the data in Fig. prefers tot the drag reduction as a 
function of the concentration of the rod-like polymers, 
presented in Prandtl-Karman coordinates. The results 
of the theoretical predictions are shown in Fig. 1161 and 
the qualitative agreement with the experimental obser- 
vations is obvious. Another interesting comparison with 
exp erimental findings is availa ble due the work presented 
in ( Amarouchene et all l2007f ) where the percentage of 
drag enhancement and reduction were measure as a func- 
tion of v p . The quantitative comparison needs a careful 
identification of the material parameters in the theory 
and the experiment, and th is was described in detail in 
( Amarouchene et all l2007h . Fig. [17] shows the com- 
parison of the percentage of drag reduction (enhance- 
ment) between the theoretical predictions and the ex- 
perimental results. The two data sets shown pertain to 
c.p = 250wppm and c p = 500wppm. The agreement be- 
tween theory and experiment is very satisfactory. 



VqS = p'L 



(6.7) 



VII. DRAG REDUCTION BY BUBBLES 



Usually in experiments the system size and the working 
fluid remain the same. Therefore, vq and L are constants 



In order not to leave the reader with the impres- 
sion that polymers are the only additives that can re- 
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FIG. 16 f~ 1/2 as the function of log 10 (7?e/ 1/2 ) with A = 1 
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of minute bubbles is relatively straightforward, since un- 
der such conditions the bubbles only renormalizc the den- 
sity and the viscosity of the fluid, and a one-flu i d mode l 
suffices to describe the dynamics (jL'vov et al ], l2005brl . 
The fluid remains incompressible, and the equations of 
motion are basically the same as for a Newtonian fluid 
with renormalized properties. The amount of drag re- 
duction under such conditions is however limited. But 
when the bubbles increase in size, the one-fluid model 
loses its validity since the bubbles become dynamical in 
the sense that they are no longer Lagrangian particles, 
their velocity is no longer the fluid velocity at their cen- 
ter, and they begin to fluctuate under the influence of 
local pressure variations. The fluctuations of the bubbles 
are of two types: 1) the bubbles are no longer spheri- 
cal, distorting their shape according to the pressure vari- 
ations, and 2) the bubbles can oscillate radially (keep- 
ing their spherical shape) due to the compressibility of 
the gas inside the bubble. The first effect was stud- 
ied numerically using th e "fron t tracking" algor ithm of 
iKawamura and Kodamal (|2002h : iLu et all (|2005l ). How- 
ever, the results indicate either a drag enhancement, or 
a limited and transient drag reduction. This leads one to 
study the possibility of explaining bubbly drag reduction 
by bubbl e oscilla t ions. Indeed, a theoretical model pro- 
posed by iLegnerl (1984) successfully explained the bub- 
bly drag reduction by relating turbulent viscosity in the 
bubbly flow to the bulk viscosity of the bubbles. While 
the bulk viscosity is important only when the bubbles 
are compressible, it is important and interesting to see 
how and why it affects the charactistics of the flow. Here 
we review how drag reduction is caused by bubbles when 
bubble oscillations are dom inant. F i nally we compare our 
finding with the results of iLegne r (1984), showing that 
a nonphysical aspect of that theory is removed, while a 
good agreement with experiment is retained. 
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FIG. 17 Percentage drag reduction as a function of Re for 
two values of the concentration of the rod-like polymer. Sym- 
bols are results of experi ments and solid lines-theory, cf. 
(jAmarouchene et al.1 . 120071 ) 



duce the drag, or that they provide the only techno- 
logically preferred method, we discuss briefly drag re- 
duction by other addi t ives l ike surfactants and bubbles 
( Gvr and Bewersdorfll |l995). Generally speaking, the 
understanding of drag reduction by these additives lags 
behind what had been achieved for polymers. The im- 
portance of drag reduction by bubbles cannot be however 
overestimated; for practical applications in the shipping 
industry the use of polymers is out of the question for 
economic and environmental reasons, but air bubbles are 



potentially very attractive ( Kitagawa et all l2005f) . 
The theory of drag reduction by small concentrations 



A. Average equations for bubbly flows: the additional 
stress tensor 

Dealing with bubbles we cannot put the fluid den- 
sity to unity any longer and we must display it explic- 
itly. Here a Newtonian fluid with density p is laden 
with bubbles of density p B , and radius R which is much 
smaller than the outer scale of turbulence L. The vol- 
ume fraction of bubbles C is taken sufficiently small 
such that the direct interactions between bubbles can 
be neglected. In writing the governing equations for 
the bubbly flow we will assume that the length scales 
of interest are larger than the bubble radius. Later 
we will distinguish however between the case of mi- 
crobubbles whose radius is smaller than the Kolmogorov 
scale r] and bubbles whose radius is of the order of rj 
or slightly larger. For length scales larger than the 
bubbles one writes dBiesheuvel and Wiingaarden . 1984 : 



ISanghai and Didwania . 1993t IZhang and Prosperetti 
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1994): 



r dU 



p{l-C)[—+U-VU] = -Vp+V-a + CV-T . (7.1) 



d(l - C) 
di 



V-(l-C)E/ = 



(7.2) 



In these equations, U is the velocity of the carrier fluid, 
and 



p ( dUi w 3 



(7.3) 



Here po = pvo- The effect of the bubbles appears in two 
ways, one through the normalization of the fluid density, 
p — > p(l— C), and the second is via the extra stress tensor 
r. Note that in the case of polymers we accounted only 
for the additional stress tensor since the actual volume 
fraction of polymers is minute, typically less then 10~ 3 . 
A typical value of the bubbles volume fraction is 0.05 or 
even more, and cannot be neglected. 

The evaluation of the extra stress tensor r is tedious. 
There are a number of forces acting on bubbles, but fol- 
lowing lLo et al. I (120061) we neglect gravity and the lift 
force. In ( Lo et alT 20061 ) it was found that the add- 
mass force due to bubble oscillations and the add-mass 
force due to bubble acceleration are of the same order of 
magnitude. The viscous forces were of course taken into 
account. The result is of the following form: 



L w - U) ■ (w - U) - RR - | R 



y2 



~(w-U)(w-U) + ^fM,a- 



}■ 



(7.4) 



Here J is the unit tensor, w is the bubble velocity which 
differs from the carrier velocity U. R stands for the rate 
of change of the bubble radius R. The relative impor- 
tance of the various terms in r depends on the values 
of Re and also on the so called Weber number, which is 
defined as 



We 



p\w- U\ 2 R 

7 



(7.5) 



The Weber number is the ratio of the kinetic energy of the 
fluid associated with the bubble motion over the surface 
energy of the bubble itself due to surface tension 7. When 
We<§; 1 the bubbles can be considered as rigid spheres. If 
We is sufficiently large, the bubbles begin to deform and 
oscillate; this contributes a significant contribution to r. 
In the following section we show that this can be crucial 
for drag reduction. 

For very small bubbles (micro-bubbles) of very small 
density the last term in t , i.e. pqS which is the viscous 
contribution, is the only one that survives. When this is 
the case the bubble contribution to the stress tensor can 
be combined with <x in Eq, (|7.ip , resulting in the effective 
viscosity given by 



The study of drag reduction under this renormaliza- 
tion of the viscosity and the density was presented in 
( L'vov et al ., 2005b) with the result that drag reduction 
can be obtained by putting the bubbles out of the viscous 
sub-layer and not too far from the wall. The amount of 
drag reduction is however rather limited in such circum- 
stances. 



B. Balance equations in the turbulent boundary layer 

At this point we apply the formalism detailed above to 
the question of drag reduction by bubbles in a stationary 
turbulent boundary layer with plain geometry. This can 
be a pressure driven turbulent channel flow or a plain 
Couette flow , which is close to the circ ular Couette flow 
realized in ( van den Berg et all 120051 ). As before , we 
take the smallest geometric scale to be 2L, the unit vector 
in the streamwise and spanwise directions be x and z 
respectively, and the distance to the nearest wall be y <JC 
L. The variables that enter the analysis are modified by 
the density, except for the mean shear, which remains 
as in Eq. (|2.2[) . The turbulent kinetic energy density is 
modified to read 



K, = \p{l-C){\u\*) 



and the Reynolds stress 



W=-p{l-C)(U X Uy) 



(7.7) 



(7. 



1. Momentum balance 



From Eq. (|7.ip we derive the exact equation for the 
momentum balance by averaging and integrating in the 
usual way, and find for i/«L: 



pP = p S + W + C(T xy ) 



(7.9) 



Here P is the momentum flux toward the wall. 

For C = Eq. (O is the usual Eq. ([3~T3|) satisfied 
by Newtonian fluids. To expose the consequences of the 
bubbles we notice that the diagonal part of the bubble 
stress tensor r [the first line in the RHS of Eq. (|7.4[) ] 
does not contribute to Eq. (|7.9[) . The xy component of 
the off-diagonal part of r is given by the 2nd line in Eq. 
(|7.4p . Define the dimensionless ratio 



c 



(K - U X )(Wy ~ Uy)) 
2(U X Uy) 



(7.10) 



For later purposes it is important to as sess the s i ze an d 
sign of C- This object was analyzed in ( Lo et all . l2006h . 
and it was shown that for small values of Re, £ is small. 
On the other hand, for large Re the fluctuating part of 
w is closely related to the fluctuating part of it. The 
relation is 



PeS 



(l+oC)p 



(7.6) 



(7.11) 
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This implies that (k2 and is positive definite, as we in- 
deed assume bellow. With this definition we can simplify 
the appearance of Eq. (|7.9p : 



(7.12) 



1 + C(C-1) 1A , 

= H W , 



1 -C 



with /Lt e g defined by (|7.6 



2. Energy balance 

Next, we consider the balance of turbulent energy in 
the log-layer. In this region, the production and dissipa- 
tion of turbulent kinetic energy is almost balanced. The 
production can be calculated exactly, WS. The dissipa- 
tion of the turbulent energy is modeled by the energy 
flux which is the kinetic energy IC(y) divided by the typ- 
ical eddy turn over time at a distance y from the wall, 
which is y/ p(l — (Tjy/bVfc where 6 is a dimensionless 
number of the order of unity. Thus the flux is written 
as bK?!" 2 1 yyj p(\ — C). The extra dissipation due to the 
bubble is C(TijSij) where = dui/dxj. In summary, 
the turbulent energy balance equation is then written as: 



bK?l 2 



V7>(1 " C)y 



+ C(njSij) = WS 



(7.13) 



As usual, the energy and momentum balance equations 
do not close the problem, and we need an additional re- 
lation between the objects of the theory, and we define 
c B in a manner following Eq. (|2.6p : 



W = c 2 K 



(7.14) 



Clearly, ^lim c B = c N and for small C (noninteracting 

bubbles) c B — cx C . Having in mind that in numerical 
simulat ions of the incompressible bubbly flow, it w as re- 
ported (jBellakhel et. al.l . 120041 : lLa 



^ance et. al~l . ll99lh that 



c B is slightly smaller than its Newtonian counterpart we 
can write 



(7.15) 



with positive coefficient (3 of the order of unity. We are 
not aware of direct measurements of this form in bubbly 
flows, but it appears natural to assume that the param- 
eter (3 is y-independent in the turbulent log-law region. 



C. Drag reduction in bubbly flows 

In this section we argue that bubble oscillations are 
crucial in enhancing the effect of drag reduction. This 
con clusion is in line w i th the experimental observation 
of ( van den Berg et all 120051 ) where bubbles and glass 
spheres were used under similar experimental conditions. 
Evidently, bubble deformations can lead to the compress- 
ibility of the bubbly mixture. This is in agreement with 



the simulation of ( Ferrante and Elgobashil |2004 ) where 
a strong correlation between compressibility and drag re- 
duction were found. 

To make the point clear we start with the analysis of 
the energy balance equation (]T. 13[) . The additional stress 
tensor Eq. (|7.4p has a diagonal and an off diagonal 
part. The off-diagonal part has a viscous part that is 
negligible for high Re. The other term can be evaluated 
using the estimate (|7.11l) . leading to a contribution to 
(r,-,'S«) written as 



{-(w-U)(w-U) : Vu) 



2(uu : Vtt) 



(7.16) 



The expression on the RHS is nothing but the spatial 
turbulent energy flux which is known to be very small 
in the log-layer compared to the production term on the 
RHS of Eq. (f7TT3|) . We will therefore neglect the off- 
diagonal part of the stress tensor in the energy equation. 
The analysis of the diagonal part of the stress depends on 
the issue of bubble oscillations and we therefore discuss 
separately oscillating bubbles and rigid spheres. 



1. Drag reduction with rigid spheres 

Consider first situations in which R — 0. This is the 
case for bubbles at small We, or when the bubbles are 
replaced by s ome particles which are le ss dense than the 
carrier fluid ( van den Berg et all l2005h . When the vol- 
ume of the bubbles is fixed, the incompressibility condi- 
tion for the Newtonian fluid is unchanged, and su = 0. 
The diagonal part of t, Eq. (|7.4p . due to the incompress- 



ibility condition sa — 0, has no contribution to (TySy). 
Due to the discussion after Eq. (|7.16p the term due to the 
extra stress tensor can be neglected and the energy bal- 
ance equation is then unchanged compared to the New- 
tonian fluid. Note that the momentum balance equation 
is nevertheless affected by the bubbles. Putting (|7.14p 
into (|7.13p . we have after simple algebra 



S 2 y 2 c 6 



W = p(l-C) 



(7.17) 



To assess the a mount of drag reducti on we will consider 
an experiment (jKitagawa et al.. 2005) in which the veloc- 
ity profile (and thus S) is maintained fixed. Drag reduc- 
tion is then measured by the reduction in the momentum 
flux P. Substituting Eq. ([7TT71) into Eq. ([7TT2"|) in which 
one can neglect p e gS, and replacing S(y)y = 1/k k we 
find 



P = 



P (l-C + CC)cl 



:b 2 



(7.18) 



If there are no bubbles (C 
turn flux P., reads 



0), the Newtonian momen- 



[h 2 



(7.19) 



The percentage of drag reduction can be defined as 



%DR = 



P =1 (l-C + (C)cl 



l-C + 3/3 )C 



(7.20) 



At small Re, ( = and the amount of drag reduction 
increases linearly with C. If Re is very large we ex- 
pect £ ss 2, and then the drag is enhanced. This result 
is in pleasing agreement w ith the experimental data in 
van den Berg et"aill2005h . Indeed, the addition of glass 



beads with density less than water caused drag reduction 
when Re is small, whereas at Re ~ (10 6 ), the drag was 
slightly enhanced. 



2. Drag reduction with flexible bubbles 

If the value of We is sufficiently large such that R ^ 0, 
the velocity field is no longer divergenceless. Also, to 
evaluate the extra stress tensor one needs the explicit 
dynamical equat ion for the bubble radiu s . Thi s equation 
was provided by IZhang and Prosperettil ( 1994) . and was 
analyzed carefully in (|Lo et all l2006h . Using the results 
of the analysis in Eq. (|7.4[) one can simplify the extra 
stress tensor (for calculating the correlation (lyS^))) to 
the form 



P 



4 1 

-/iV • ul + -{w - U)(w - U) 



(7.21) 



The extra turbulent dissipation due to the bubble is 
(rijSij). In light of the smallness of the term in Eq. (|7.16[) 
we find 



' y °y / 



= ^ 



(7.22) 



The term |/xs 2 j is of the same form as the usual dissipa- 
tion term /xSy-Sy and therefore we write this as: 



t\c>)-A p W 3/2 -A V 



(7.23) 



where A is an empirical constant. Finally, the energy 
equation becomes 



b(l-C)+AC K?l 2 



pVT^C 



y 



ws 



(7.24) 



As we derived Eq. (|7.19[) . we again specialize the situa- 
tion to an experiment in which S is constant, and com- 
pute the momentum flux 



P = 



(i-c) 2 (i-c + cc) 4 



(1-C + fC) 2 4& 2 
The degree of drag reduction is then 

^2 



%DR = 1 - 



(1-C) 2 (l-C + (C) c 6 



(7.25) 



(7.26) 
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FIG. 18 Predicted values of drag reduction with a = 2 and 
different values of A. In a dashed line we reproduce the predic- 
tions of Legner's theory which suffer from an unphysical drag 
enhancement at C — 0, For A = (rigid spheres) we find only 
increasing drag enhancement as a function of C. For small 
values of A we have first a slight drag enhancement, and then 
modest drag reduction. For large values of A, associated with 
strong bubble oscillations, we find significant values of drag 
reduction. 



Note that A is an unknown parameter that should de- 
pend on We, and so its value is different in different ex- 
periments. The percentage of drag reduction for various 
values of A are shown in Fig. [18] where we chose a = 2 
and for simplicity we estimate c B = c N . One sees that for 
a = 2 and ^4 = (where the latter is associated with rigid 
bubbles), we only find drag enhancement. For small value 
of A, or small amplitudes of oscillations, small concentra- 
tions of bubbles lead (for a = 2) to drag enhancement, 
but upon increasing the concentration we find modest 
drag reduction. Larger values of A lead to considerably 
large degrees of drag reduction. For A — 0.15, the re- 
sult agrees reasonably with L egner's theory which pre- 
dicts %DR w 1- 5(1- C) 2 /4 (|Legnerl . ll98l . Note that 
according to Legner, there should be considerable drag 
enhancement when C = 0. This is of course a nonsen- 
sical result that is absent in our theory. For A = 0.8, 
%DR « 4C for small C. This is the best fit to the exper- 
imen tal results which are reported in ( Kitagawa et all 
120051) . 



VIII. SUMMARY AND DISCUSSION 

We offered a review of drag reduction by polymers and 
bubbles. Quite generally, we have shown that drag reduc- 
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tion can be understood in great detail by using only a few 
equations that govern the budgets of energy and momen- 
tum. Both polymers and bubbles open up an additional 
channel for dissipation, and thus pose a fundamental rid- 
dle: why do they reduce the drag. The answer in all 
cases is fundamentally the same: the same agents reduce 
the momentum flux from the bulk to the wall, and this 
is the main effect leading to drag reduction. The re- 
duction in the momentum flux overwhelms the increase 
dissipation. We stressed above, and we re-iterate here, 
that this mechanism depends on the existence of a wall 
which breaks the translational symmetry. Drag reduction 
must be discussed in the context of wall-bounded flows to 
make full sense. Of course, additives may influence also 
the spectrum of turbulent fluctuations in homogeneous 
and isotropic turbulence, but this is another story, quite 
independent of drag reduction. 

One should make a great distinction between the two 
cases. The phenomenon of drag reduction by polymers 
exhibits interesting universal features which are shared 
even by flexible and rod-like polymers. These features are 
the most prominent experimental results that required 
theoretical understanding. We have shown that the MDR 
has a special significance in being an edge solution of tur- 
bulent flows. Trying to increase the reduction of the drag 
behind what is afforded by the MDR would re-laminarize 
the flow. This may be the central theoretical insight that 
is offered in this review, providing a simple and intuitive 
meaning to the nature of the MDR. This explains why 
flexible and rod-like polymers have the same MDR, even 
though they approach the MDR in distinctly different 
ways. 

Once the theory was put forth to explain the universal 
aspects of drag reduction by polymers, it became also 
clear that it can easily explain, in considerable quanti- 
tative detail, also the non-universal aspects, including 
cross-overs due to small concentrations of polymers, low 
values of Tie, and small number of monomers in the poly- 
mer chains. We trust that this detailed understanding 
can help in designing and optimizing the use of polymers 
in practical applications of drag reduction. 

On the other hand, the case of bubbles exhibits much 
less universality, since the placement of the bubbles with 
respect to the wall and their actual density profile have 
crucial consequences regarding their efficacy as drag re- 
ducers. The main conclusion of our study is that bubble 
oscillations can contribute decisively to drag reduction 
in turbulen t flows. In agreement wit h the experimental 
findings of ( van den Berg et al.[ |2005r). we find that rigid 
bubbles tend to drag enhance, and the introduction of os- 
cillations whose amplitude is measured by the parameter 
A (Fig. [T8|) increases the efficacy of drag reduction. The 
main drawback of the present study is that the bubble 
concentration was taken uniform in the flow. In real- 
ity a profile of bubble concentration may lead to even 
stronger drag reduction if placed correctly with respect 
to the wall. A consistent study of this possibility calls for 
the consideration of buoyancy and the self-consistent so- 



lution of the bubble concentration profile. Such an effort 
is beyond the scope of this review and must await future 
progress. 
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APPENDIX A: The Hidden Symmetry of the Balance 
Equations 

Consider the following identity: 

v + {y + ) = 1 + a(y+ - 5+) 

= [1 + a(y+ - S) + a{6 - 6+)} 
a 



1 + 



where 



g(8) = 1 + a(5 - 5+) , 5 > 5+ 



(Al) 



(A2) 



Next introduce newly renormalized units using the effec- 
tive viscosity g(6), i.e. 



V ~ 9(5)' -g(5) 



S* = S+g{6) , W t = W + . 

In terms of these variables the balance equations are 
rewritten as 



[1 + a(y* - + W* = 1 



[1 + a(yi - 



A» 



= s* 



(A4) 
(A5) 



These equations are isomorphic to (13.19)) and (|3.20p with 
5 + replaced by 6*. The ansatz (|3.21[) is then replaced 
by A(a) = S + g(S)~ 1 f(aS*) This form is dictated by the 
following considerations: (i) A (a) — > S + when a — ► 0, 
(ii) all lengths scales in the re-scaled units are divided 
by g{S), and thus the pre- factor in front of / becomes 
5 + /g(5), and (iii) ad + in Eq. (|3.18|) is now replaced in 
Eq. (|A5|) by a$i , leading to the new argument of /. Since 
the function A (a) cannot change due to the change of 
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variables, the function A (a) should be identical to that 
given by Eq. (|3~2"Tj) : 

S + f(a6+) = • (A6) 

Using the explicit form of g(S) Eq. (|A2|) . and choosing 
(formally first) 6 = 6* = we find that /(£) = 1/(1 - £). 
It is easy to verify that this is indeed the solution of 
the above equation for any value of 5* , and therefore the 
unique form of Eq. (|3.22[) . 
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